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ABSTRACT 

The GUI red giant star e Oph has been found to exhibit several modes of oscillation by the MOST mission. We interpret the observed 
frequencies of oscillation in terms of theoretical radial p-mode frequencies of stellar models. Evolutionary models of this star, in both 
shell H-burning and core He-burning phases of evolution, are constructed using as constraints a combination of measurements from 
classical ground-based observations (for luminosity, temperature, and chemical composition) and seismic observations from MOST. 
Radial frequencies of models in either evolutionary phase can reproduce the observed frequency spectrum of e Oph almost equally 
well. The best-fit models indicate a mass in the range of 1.85 ± 0.05M o with radius of 10.55 ± 0. 15/? Q . We also obtain an independent 
estimate of the radius of e Oph using high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR 
instrument. The measured limb darkened disk angular diameter of e Oph is 2.961 ± 0.007 mas. Together with the Hipparcos parallax, 
this translates into a photospheric radius of R = 10.39 ± 0.07R o . The radius obtained from the asteroseismic analysis matches the 
interferometric value quite closely even though the radius was not constrained during the modelling. 
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1. Introduction 



Stars: oscillations - Stars: interiors - Stars: fundamental parameters - Techniques: 



Asteroseismology of red giant stars has, in recent years, taken 
a leap forward with the discovery of pulsation s in several G- 
and K-type giant stars, bo th from the ground (Frandse n et al.l 
2002; de Ridd er et al.ll2006l) and from space dBarban et al.l 2007; 
Kall inger et al.ll2008bt iHekker et al]|2008l) . Oscillations in the 
G giant e Oph (HD 146791, HR 6075, HIP 79882) was 
first detected in spect roscopic observations from the ground 
dde Ridder et al.l 120061) . although the average large separation 
could not be distinguished between two possible values due to 
the daily alias problem. Sub sequent observations by the MOST 
satellite (Wal ker et al.l 12003) led to the discovery of at least 9 
radial modes with a n average large separation of 5.3 +0.1 fiHz 
dBarban et al.ll2007l) . 

This work makes an attempt to interpret the observed fre- 
quencies of e Oph in terms of adiabatic oscillation modes of stel- 
lar models in the relevant part of the HR diagram. We construct 
red giant models in which the luminosity is provided by either 
hydrogen burning in a shell outside the helium core, or both shell 
hydrogen burning and core helium burning. We make quantita- 
tive comparisons of these models to the MOST frequencies of 
e Oph to determine the stellar parameters like mass, age, radius, 
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and chemical composition. iKallinger et alj d2008al) have earlier 
presented stellar models in the shell hydrogen-burning phase for 
e Oph, based on asteroseismic data. A similar study of oscilla- 
tions in the red gia nt £ Hya in terms of helium b urning models 
was carried out bv lChristensen-Dalsgaardl d2004l) . 

Interferometric measurements of stellar radii are particu- 
larly discriminating for models, in particular when combined 
with asterose i smic frequ encies, as noticed, for instance, by 
ICreevev et ail d2007l) and ICunha et alj d2007l) . In this work we 
report a new interferometric determination of the radius of 
e Oph. While the direct measurement of the radius of a red giant 
is a useful result in itself, for e Oph it provides the first oppor- 
tunity of testing the relevance of theoretical models for red gi- 
ants which have been calibrated with asteroseismic input. In this 
study, neither did we use the interferometric radius as an input 
to the modelling, nor did the interferometric analysis draw upon 
the asteroseismic information in any way. Thus the radius of the 
stellar models that fit the seismic data best can be tested against 
the independently measured interferometric radius. 

In Sect.|2]we describe the details of the stellar models that we 
have constructed and in Sect. [3] we compare the theoretical fre- 
quencies obtained from these models with the observed MOST 
frequencies. In Sect. [4] we present our new interferometric mea- 
surement of the angular diameter of e Oph. In Sect. [5] we com- 
pare the radii of our best seismic models with the interferometric 
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measurement, and discuss our results with similar studies carried 
out earlier. 



2. Stellar models and theoretical frequencies 

We constructed a grid of stellar models with various input pa- 
rameters using the Y ale Rotating Evolutionary Code (YREC) 
dGuenther et al.llT992b . This code is capable of producing con- 
sistent stellar models for low mass giant stars bot h in the shell H- 
burning phase and the core He-burning phase (Demarqu e et al.l 
2008). We describe these two sets of models below. The radial 
and nonradial pulsation freq uencies of each model are calculated 
by the oscillation code JIG (iGuenfherl 19941) . 



2.1. Input physics 



The models use the late st 
( [Rogers & Navfonoyl |2002|) 
(Iglesias & Rogers! 1996|). supplemen te d by the 



OPAL 
and 



equation of state 
OPAL opacities 
low tem- 



perature opacities of Ferguson et al. (2005). The nuclear 



reaction rates from Bahcall & Pinsonneault (1992) are used 



Diffusion of helium and heavy elements have been ignored in 
the post-main sequence phase of evolution. This is reasonable, 
since dredge-up by the deep convective envelope present in 
red giants would mask any effect of diffusion of elements in 
early phases. The current treatment of convectio n in the models 
is thr ough the standard mixing length theory (Bohm-Vitense 
Il958l) . which does not properly include the effects of turbu- 
lence in the outer layers, and this might substantially affec t 
the frequencies of oscillation (see e.g., IStraka et alj [2007). 
Fortunately, the uncertainty induced in the large separations 
is much less than in the actual frequencies themselves. Mass 
loss on the giant branch was not included in the calculations. 
Most of the mass loss is believed to take place quiescently as 
the star approaches the t ip of the giant branch, as is th e case in 
the c ommonly adopted Rei mersl d 19771) formulation dYi et alj 
120031) . Such mass loss does not affect the thermodynamics of 
the deep interior appreciably. Most importantly, it takes place at 
high luminosities, beyond the luminosity of e Oph on the giant 
branch. The neutrino lo sses in the core were taken from the 
work of lltoh et al.l d!989t) . 



2.2. Range of parameters 

The range of input parameters chosen for the models is dic- 
tated by the position of e Oph on the HR diagram and its 
estimated chemical composition. The values of effective tem- 
perature (log r e ff = [3.680,3.698]), luminosity (logL/L = 
[1.732, 1.8061 ), and metallicity gFe /Hl = [-0.07, -0.17]) are 
adopted from lde Ridder et al.l d2006l) . who have already carried 
out a detailed survey of these parameters from the literature. 
Tracks were constructed for different (Yq, Zo) combinations with 
Yq ranging from 0.255 to 0.280, and Zo ranging from 0.012 to 
0.015, corresponding to [Fe/H] = (-0.07, -0.17). W e adopted 
the solar abundances as given by iGrevesse & Sauval d 1998b in 
converting [Fe/H] to Z values. 

For each of these combinations, two values of the mixing 
length parameter a (the ratio of mixing length to local pressure 
scale height) have been used: 1.6 and 1. 8. Note that using a sim- 
ilar version of YREC and input physics. iKalhnger et al.l d2008a) 
used a = 1 .74 in their stu dy of e Oph. This is the value of a 
adopted by lYi et al.l d2003l) for the standard solar model cali- 



bration. The radii of red giant models depend sensitively on the 
choice of a. 

For each (Yo,Zo,a), the mass has been varied between 1.8 and 
2.4M Q to check for overlap in the box. In most models, overshoot 
at the edge of the convective core present on the main sequence 
was assumed to have negligible effect on the advanced evolution. 
A few models were constructed with overshoot of 0.2 times the 
pressure scale height at the convective core edge. Because the 
size of the convective core is not too large in this mass range, the 
overall structural effect of core overshoot is modest; but evolu- 
tionary timescales are slightly increased w hen core overshoot is 
taken into account dDemarque et al.ll2004l) . 



2.3. Shell H -burning models 

Our first set of models for e Oph are on the ascending red giant 
branch. These models are characterised by an inert helium core 
surrounded by a thin hydrogen burning shell. The mass of the 
shell varies between 0.012M o and 0.0007 M depending on the 
mass and age. The mass of the shell decreases as the star ascends 
the red giant branch. We concentrate on models that lie within 
the e Oph box on the HR diagram. For each evolutionary track 
that traverses the box, several models at slightly different ages 
are constructed so as to span the box. The theoretical frequencies 
of these models are compared to the observed frequencies of 
e Oph in Sect. |3] 

Each evolutionary track was started in the pre-main sequence 
phase and evolved continuously through core hydrogen burning 
and eventually shell hydrogen burning along the giant branch. 
Some of the tracks are plotted in the HR diagram in Fig. [TJ to- 
gether with the e Oph error box. The size of this box is such 
that the range of mass of models with a given set of {Yq, Zo) and 
a values that pass through the box is ~ 0.5M o (see top panel 
of Fig. [T}. Typically, the mass lies between 1.8M Q and 2.4M Q , 
depending on the values of the other parameters. The shift in 
the tracks with these parameters is also significant (see bottom 
panel of Fig. [T}. Since the outer convective layer of stars in this 
mass range is extremely thin during the main sequence phase, 
the tracks with different values of a are almost identical in that 
phase (e.g., tracks C and D in Fig.[TJ. But the shift of the track on 
the giant branch is significant because of the extended convec- 
tive envelope. In most cases in this study, however, frequencies 
of models inside the e Oph box with a = 1.6 had a poor match 
with the observed frequencies. Since the tracks move redwards 
with decreasing a, this can be traced to larger mass (and hence 
larger radii inside the box) of the a = 1.6 models compared 
to the a — 1.8 models. Keeping in mind the adopted range of 
[Fe/H] for the star, the metallicity of the initial model can only 
be varied between 0.012 and 0.015 with corresponding appropri- 
ate change of initial helium abundance between 0.255 and 0.280 
to span the entire e Oph box on the HR diagram. 



2.4. Core He-burning models 

The second set of models are in a later stage of evolution than 
the first. These models have helium burning in the core of the 
star. Hydrogen burning in a thin shell outside the core is also 
present. Considering the adopted metallicity value for e Oph, 
these models represent the so-called "red clump" stars, rather 
than metal-poor horizontal branch stars. Indeed, the evolutionary 
tracks of the models that we constructed lie very close to the 
red giant branch, and therefore, overlap with the error box of 
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Fig. 1. Stellar evolutionary tracks from the ZAMS to the shell 
H-burning RGB are shown with respect to the position of e Oph 
on the HR diagram. The top panel shows tracks for different 
masses (M/M e indicated) with identical initial chemical com- 
position (To = 0.270, Zo = 0.012) and mixing length (a = 1.8). 
The bottom panel shows tracks with the same mass (M = 2.2M ) 
but with different values of Yq, Zo, and a. 



e Oph on the HR diagram. Fig. [2] illustrates these tracks in the 
HR diagram. 

For the range of chemical composition used in our models, 
it turns out that the helium ignition in the core at the tip of the 
red giant branch takes place under degenerate conditions for the 
lower mass stars (M/M Q $ 2.1). This is the well-known "he- 
lium fl ash" mechanism first studied by [Schwarzschild & H arm 
(fl962h . 

For stars of slightly higher mass (M/M e ^ 2.1), however, 
the core is not degenerate at the instant helium burning temper- 
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Fig. 2. Stellar evolutionary tracks from the ZAHB to the end of 
He-burning main sequence are shown with respect to the posi- 
tion of e Oph on the HR diagram. The direction of evolution is 
from the lower left to the upper right corner of the graphs. The 
tracks in each panel are the successors of those in Fig. Q] The 
top panel shows tracks for different masses (M/M indicated) 
with identical values of initial chemical composition and mix- 
ing length. The bottom panel shows tracks with the same mass 
(M = 2.2M ) but with different values of To, Zo, and a (given in 
Fig.0. 



atures are reached, and therefore, helium ignition takes place in 
a controlled fashion. In this latter case, it is numerically easy to 
continue the evolution of the model past the tip of the red giant 
branch, and onto the red clump phase. 

For the helium flash scenario, however, the numerical stabil- 
ity of the evolution code at the tip of the giant branch is far less, 
and only a computationally expensive algorithm involving sub- 
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tie handling of various parameters can guide the model past the 
runaway helium ign ition process and settle it on t o a stable phase 
of hel ium burning (Demarque & Mengel [i97ll) . iPiersanti et al.l 
(2004) have demonstrated that even for stars which undergo a 
violent helium flash, the subsequent evolution of a model on 
the helium-burning main sequence, i.e., once stable core helium 
burning has been established, is not very sensitive to the prior 
history of helium ignition. Specifically, the behaviour of mod- 
els which have been evolved from appropriate zero age hori- 
zontal branch (ZAHB) models with quiescent helium burning 
in the core is remarkably similar to that of models which have 
been actually evolved through the helium flash phenomenon. Of 
course, the make-up of the ZAHB model is crucial - it must re- 
flect the properties of a model that has settled on the helium- 
burning main sequence after having gone through the helium 
flash. The critical factor in this starting model, apart from the 
total mass and the chemical composition at the tip of the giant 
branch , is the mass of the heli um core. Hydrodynamical st udies 
in 2D (ICole & Deupredl 1983b and 3D (iMocak et al.ll2008l) con- 
firm this picture, except for possible second order mixing effects 
due to turbulent overshoot at the convective-radiative interface, 
which cannot, at this point, be estimated precisely. 

For the low mass models we have, therefore, followed the 
approach of circumventing the numerical difficulties encoun- 
tered in handling the helium flash, as done b y most authors 
(e.g., iLee & Demarque Ill990t ISweigartlll987l) . For each mass 
and chemical composition evolution was continued on the red 
giant branch till the onset of helium flash. The evolution of the 
red clump model was then re-started from a ZAHB model of the 
same total mass with identical chemical composition and helium 
core mass as the corresponding model at the onset of helium 
flash at the tip of the red giant branch. 

Fig. H] illustrates the behaviour of the central regions of a 
2M Q model near the tip of the giant branch, where helium igni- 
tion takes place. The four panels describe the change, as a func- 
tion of time, or equivalently, as a function of maximum tem- 
perature in the star, T m - dx , of the following quantities: the mass 
contained interior to the shell at T max , w(7 , max )/M , the degen- 
eracy parameter 77, the energy generation rate due to the triple- 
alpha reaction, e(3a), and finally the mass of the helium core, 
'"He-core /M Q . Note that the shell with the highest temperature is 
not central (i.e., m(T m . dx )/M Q + 0) and changes with time. This 
is because of neutrino losses, which are most effective at the 
higher densities near the centre. As the degeneracy increases, 
neutrino cooling causes such an inversion of the temperature 
profile near the centre till the helium ignition temperature is 
reached. Such off-centre ignition of helium is typical in degen- 
erate helium cores. The energy released in the helium ignition 
reduces the degeneracy, and the shell of maximum temperature 
moves back to the centre of the star. The degeneracy parameter 
77 is a measure of the degree of degeneracy of the electron gas. 
It is a dimensionless parameter used to quantify the relationship 
between the electron density and pressure in a partially degen- 
erate electron ga s. The detailed f ormalism used in the models is 
that described by IClavtonl (1 1 9681 p. 64). We note, while inspect- 
ing Fig. 3, that the quantity r\ varies between -00 in the ideal gas 
case and +00 in a fully degenerate Fermi-Dirac gas. The mass of 
the helium core keeps increasing throughout the red giant phase 
due to hydrogen burning in the shell immediately above it till the 
onset of helium burning. The maximum mass of a degenerate he- 
lium core at helium flash is typically ~ 0.45M o , irrespective of 
the total mass of the star, and depends slightly on the other stellar 
parameters. 
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log T max (in K) 



7.95 



Fig. 3. The change, as a function of time (upper x-axis), or equiv- 
alently as a function of maximum temperature in the star, r max 
(lower x-axis), of (a) the mass interior to the shell at T max , 
m(J mdx )IM Q , (b) the degeneracy parameter, 77, (c) the energy 
generation due to the triple-alpha reaction, e(3ch), and (d) the 
mass of the helium core, TWee-core/Ms, close to the onset of He- 
burning of a 2M Q star are shown. The red solid line denotes the 
pre-He-ignition phase, while the blue dotted line shows the post- 
He-ignition phase. 



The actual age of the model on the helium-burning main 
sequence cannot, of course, be assigned accurately due to the 
"missing" period of the helium flash. However, the total dura- 
tion of the helium flash phenomenon and the subsequent stabil- 
isation of t he star on the he lium-burning main sequence is only 
~ 1.5 Myr (IPiersanti et al.l 120041) . and hence the uncertainty in 
the age in our helium-burning models is quite small. 

For slightly higher mass models (M/M G 2.2), the evo- 
lution of the star is followed continuously from the pre-main 
sequence stage till the red clump stage. Since the helium igni- 
tion at the tip of the giant branch occurs under non-degenerate 
conditions, there are no numerical problems in such cases. The 
transition mass from violent to quiescent helium ignition is a 
fu nction of chemical co mposition. It has been studied in detail 
bv lSweieart et al.J dl989). More recent illustration is found in the 
tracks of lYi et all (120031). which were co nstructed with a similar 
version of YREC (lDemarque 'et~ai]|2008l) . 
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Fig. 4. The different timescales of the stellar evolutionary tracks 
in each crossing of the e Oph box for the shell H-burning, post 
He-ignition, and core He-burning main sequence are shown. In 
the top panel the full evolutionary track of a 2.3M Q model is 
shown on the HR diagram. The bottom panel shows the lumi- 
nosity of the star as a function of its age. The luminosity limits 
of e Oph are indicated by the two horizontal lines, and the time 
spent by the star within these limits in each phase are indicated 
by the projections on the x-axis. 



There is an uncertainty in the age of the core helium-burning 
models due to the assumption of no mass loss on the red giant 
branch. Since the amount of mass lost in the giant phase is not 
known, the mass of the ZAHB model cannot be assigned accu- 
rately. Therefore, the ages of helium-burning models given in 
Tables Q] and [4] are, in fact, upper limits to the real ages. 

2.5. Timescales of evolution 

A star of a given mass might cross the e Oph errorbox on the 
HR diagram three times - once upwards during shell H-burning, 
once downwards during the stabilisation of the star just after He 
ignition, and once again upwards during the phase of stable core 
He-burning. The timescales of evolution in these three phases 
are quite different, and consequently the time spent inside the 
e Oph box differs vastly. Fig.|4]illustrates the timescales of evo- 
lution in each crossing of the box for a 2.3M star. During the 



shell H-burning phase, it spends nearly 1 .09 Myrs inside the box. 
After He ignition in the core, it spends only 0.17 Myr during the 
rapid settling towards the He-burning main sequence, and finally 
it spends 30.59 Myrs during the stable core He-burning phase. 
Similar timescales are found for other stars in this mass range. 
The time spent during the shell H-burning phase is typically ~ 20 
times shorter than that during the core He-burning phase. Thus 
the probabilities of e Oph being in the corresponding phases of 
evolution are in the same ratio. 



3. Comparison of theoretical and observed 
frequencies 

The theoretical frequencies of the stellar models were compared 
with the MOST data on e Oph. Typically, in asteroseismic mod- 
elling studies, the comparison between a stellar model and the 
observed frequency data is carried out in terms of frequency 
separations, especially the large frequency separations, rather 
than the frequencies themselves. This is done to eliminate the 
uncertainty in the theoretical absolute frequencies due to inad- 
equate modelling of the stellar surface layers. However, this is 
possible only in the happy circumstance of detection of a series 
of frequencies of the same degree and successive radial orders, 
for which the large separations can be determined. For e Oph, 
Barban et al. (2007) indeed provide the frequencies of 9 succes- 
sive radial order modes. Thus it is possible to match the observed 
large separations with the theoretical values from the models. 
However, as an additional comparison, we also match the abso- 
lute frequencies of radial modes of our stellar models with the 
MOST data. 

For each comparison, a reduced^ 2 value is computed as 



X v 



1 N 

N ^ 



VMOST - Vmodel 



1 N 

" = N- 1 S 



Avmost - Av mo del 



(JAvmost 



(1) 



(2) 



where v and Av represent the frequency and large separation re- 
spectively, and N is the number of observed modes compared 
(maximum of 9). <5vmost is the measured error in the MOST fre- 
quency and <5Avmost is the error in the large separation com- 
puted by adding the adjacent frequency errors in quadrature. 
Thus for each model, values of xl and x\ v can be computed. A 
small value of xl would usually imply a small value for^ 2 y , but 
the reverse is not necessarily true. It is possible to find models 
for which the large separations are quite close to the observed 
values, but all the frequencies are shifted by a nearly constant 
amount. Given the inadequacies of input physics in the mod- 
els, especially in the rarefied outer convective layers, the abso- 
lute values of the theoretical frequencies cannot be trusted too 
much. The large separations, on the other hand, would be rela- 
tively free from such ambiguities, and would reflect the overall 
structure of the star better. This is why we put a greater impor- 
tance on minimising %\ , rather than xl, in choosing our best 
model for e Oph. It turns out, however, that the models with 
lowes t x 2 . have reasonably s mall values for xl as well in most 
cases. iKallinger et all d20 08a) have used absolute frequencies in 
their comparison of models and observations of e Oph. 

The best-fitting models in both shell H-burning and core 
He-burning phases can reproduce the observed large separations 
fairly well. This is illustrated in Fig [5] However, it turns out that 
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Table 1. Stellar parameters for the models with lowest x\ an d;f 2 values in either shell H-burning or core He-burning phase are 
given. 
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Table 2. Calibrators used for e Oph. The limb darkened angular 
diameter #ld are given in milliarcseconds (mas). The angular 
separation y between each calibrator and e Oph is given in the 
last column, in degrees. 



50 60 
Frequency (v) (yuHz) 

Fig. 5. The large frequency separations of best-fitting radial 
modes of stellar models in the shell H-burning (red open 
squares) and core He-burning (blue filled squares) phases inside 
the e Oph errorbox are shown in comparison with the observed 
separations of e Oph. The parameters of these models are given 
in Table Q] These models are chosen on the basis of lowest x\ v 
criterion. 



the large separation corresponding to one of the observed modes 
at Vmost = 62.871 yuHz is nearly 2cr away from the theoretical 
value in all our models. This data point is always the largest con- 
tributor to;^ 2 and;^ y . We have ignored this point while choosing 
our best model. 

The parameters of the models with the lowest x\ v an d xl 
values are listed in Table [TJ Notice that the model with least x\ v 
in the shell H-burning phase has a significantly large value of 
xl due to an overall shift in the absolute frequencies. Given one 
criterion for comparison (either x\ v or xl)^ it i s clear that the 
models in the shell H-burning phase fit the data almost equally 
well as those in the core He-burning phase. Thus, the present 
data is unable to distinguish between these two phases of stellar 
evolution. However, as discussed in Sect. 12.51 the likelihood of 
the star being in the core He-burning phase is greater than it 
being in the shell H-burning phase. 

Based on the best match between observed and model large 
separations, the models indicate very similar parameters for both 
phases of evolution. We estimate the stellar parameters from not 
only the models with lowest^- 2 , but actually all models that have 
X 2 values within 50% of the least x 2 - The mass is estimated to 
be 1.85 + 0.05M o , while the metallicity of the best models are 
in the range of 0.0125 + 0.0005. The radius lies in the range of 
10.55 + 0.15/? Q . The radius of e Oph, however, can be measured 



Star 


m v 




Spect. 


6 LD (mas) 


y(°) 


HR 145085 


5.9 


2.4 


K5III 


1.677 + 0.022 


8 


HD 162468 


6.2 


3.2 


K1III 


1.154 ±0.015 


28 


HD 166460 


5.5 


2.6 


K2III 


1.439 ±0.018 


29 



independently through interferometry, as described in the next 
section. 



4. Interferometric measurements 

4. 1 . Instrumental setup 

We observed e Oph in July 2006 at the CHARA Array 
dten B rummelaar 2005) usin g FLUOR, the Fiber Linked Unit 
for Optical Recombination dCoude du Foresto et al. I I2003I) . This 
instrument is equipped with a near infrared K' band filter 
(1.9 < A < 2.35 yum). We extracted the instrumental visi- 
bilities fr om the raw data using the FLUOR data reduction 
software dCoude du Foresto et all 1 19971: iKervella et all l2004t 
Merand 2006). For all the reported observations, we used the 
CHARA baselines S2-W2, with ground lengths of 177 m, which 
is mostly a north-south baseline in orientation. The calibrator 
stars w ere chosen in the catalogue compiled by Merand et al. 
(2005), using criteria defined by these authors (Table They 
were observed immediately before or after our targets in order 
to monitor the interferometric transfer function of the instru- 
ment. For a more detailed description of the observing proce- 
d ure and the error prop agat ion, the intere sted reader is referred 
to IKervella etal J d2008l) and lPerrinl d2003l) . respectively. The re- 
sulting calibrated squared visibilities are listed in Table [3] 

4.2. Angular diameter measurement and precision 

In order to accurately measure e Oph angular diameters, we used 
many known stellar calibrators and we repeated the observation 
on two separate and consecutive nights. Since, in the end, the 
number of visibility will always be small, statistically speaking, 
the final confidence on the precision will rely more on the re- 
peatability of the the result and the consistency between the stel- 
lar calibrators. 

Our observation strategy was designed to give maximum 
precision and confidence to our results. As a result, we achieve 

- the repeatability of the result. The first night gives 0ud = 
2.881 ±0.006 mas, with areduced^- 2 of 0.9; the second night 
0ud = 2.891 + 0.005 mas, with a reduced x 2 of 0.3. This is 
consistent at the 0.01 mas level. 

- the use of multiple slightly resolved calibrators of various 
sizes. Indeed, if we assume we have an overall bias in the 
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Table 3. Squared visibility measurements obtained for e Oph are 
given. B is the projected baseline length, and "PA" is the azimuth 
of the projected baseline (counted positively from North to East). 



MJD 


S(m) 


PAD 


V 1 ± <r(V 2 ) 


53936.21641 


154.375 


-32.471 


0.03131 +0.00173 


53936.25567 


165.132 


-36.878 


0.01450 ± 0.00049 


53937.20561 


151.978 


-31.267 


0.03746 ± 0.00104 


53937.23681 


160.961 


-35.346 


0.01901 +0.00055 


53937.26478 


167.893 


-37.774 


0.01054 + 0.00028 



angular diameter estimation of the calibrators, it is going to 
lead to a differential calibrated visibility bias, depending on 
the size of the calibrator. For example, if we multiply the di- 
ameters of our calibrators by 1.05 (i.e. a 5% bias), the new 
diameter for e Oph is #ud = 2.895 + 0.005 mas, with a re- 
duced^ 2 of 2.5 instead of our result 0ud = 2.888 +0.003 mas, 
with a reduced x 2 of 1 .0. Not only our final result would 
be barely affected, moreover the reduced \ 2 becomes much 
larger, because points calibrated by our large calibrator be- 
come completely inconsistent with the rest of the batch. 

Another possible source of bias, and possibly the ultimat e one, 
is the wavelength calibration. This has been done by iMerandl 
(2005) and the process leads to a wavelength calibration ac- 
curacy of 0.005 pm. We changed the software wavelength cal- 
ibration by this amount and redid the reduction and calibration. 
This led to a bias of 0.006 mas. Hence, we added quadratically a 
0.006 mas error, leading to a final uncertainty of 0.007 mas. 



4.3. Limb darkened angular diameter and photospheric 
radius 

In order to estimate the unbiased angular diameter from the mea- 
sured visibilities it is necessary to know the intensity distribution 
of the light on the stellar disk, i.e., the limb darkening (LD). 
As we do not fully resolve e Oph, we cannot measure the limb 
darkening directl y from the data. We thus model it using the 
MARCS models dGustafsson et alJl^ Sfl for the computation 
of the intensity profile of the star, taking into account the ac- 
tual spectral tr ansmission function of the FLUOR instrument 
(Merand 2005). The LD coefficie nts have been comp uted with 
the TURBOSPECTRUM code dAlvarez & Plezlll998l) . The re- 
sult is shown in Fig. [6] It is to be noted that taking an inten- 
sity profile from a different model, say the one pr edicted by th e 
ATLAS9 model from Kurucz, using Claret's laws (IClaretll2 000), 
results in the same final result, within fractions of the statisti- 
cal error bar. The reason is that for the relatively large spectral 
bandwidth of FLUOR and considering that we measure first- 
lobe visibilities only, the difference between the MARCS and 
ATLAS models is negligible (of the order of a fraction of our 
diameter error bar). The magnitude of the limb darkening effect 
being much smaller in the infrared than in the visible, our result- 
ing limb darkened angular diameter measurement in the K band 
is largely unaffected by the choice of the limb darkening model. 

The result of the visibility fit is presented in Fig.|7]using the 
MARCS limb darkening model. We derive the following limb 
darkened disk angular diameter: 



6 hD (e Oph) = 2.961 + 0.003 (stat.) + 0.006 (syst.) mas 



(3) 



~ o.6 

0.4 
0.2 



T e( , = 5000 K, [Fe/H] = 0.0 

MARCS models: 

x log g = 2.0 
□ log g = 3.0 

--: ATLAS model, log g = 2 (Claret 2000) 



0.4 



0.2 



Fig. 6. Comparison of the intensity profiles of e Oph from the 
MARCS (re d crosses and blue squares for two values of log g) 
and Claret's (Claret 2000) (dashed curve) model. 




155 160 165 

Baseline (m) 

Fig. 7. Squared visibilities and adjusted limb darkened disk visi- 
bility model for e Oph (see also Table|3j. The lower panel shows 
the residuals to the model, with lcr (dashed line) and 3<x (dotted 
line) limits. The inset panel shows the (u, v) coverage. Note that 
the error bar given here is the formal error bar: we added some 
systematics for final result (see text). 

This value is comp atible with the spect rophotometric angular di- 
ameter estimate bv lCohen et al.l dl999l) of ft p = 3.00+0.03 mas. 
We took e Oph's parallax from the reprocessed Hipparcos cata- 
logue by Ivan Leeuwenld2007albl) : 

n(e Oph) = 30.64 + 0.20 mas (±0.65%). (4) 

This value compares well with van Alte na et al. I (fl995l) . and the 
original Hipparcos catalogue (ESAI I1997). but is more precise. 
We finally derive the photospheric linear radius: 



http://www.marcs.astro.uu.se/ 



fl(eOph) = 10.39 + O.O7fl (+0.67%). 



(5) 
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In spite of the the relatively high precision of the parallax, it is 
by far the limiting factor for the precision of the radius. 

5. Discussion 

In this study we have constructed stellar models of red giants 
in both shell H-burning and core He-burning phases and com- 
pared their adiabatic frequencies with the frequen cies of e Oph 
observ ed by the MOST satellite and published by Barba rTet al.l 
(120071) . We have also measured the radius of the star through 
optical interferometry using the CHARA/FLUOR instrument. 

We have demonstrated that the observed frequencies of 
e Oph are consistent with radial p-mode pulsations of a red gi- 
ant at the relevant position on the HR diagram. Unfortunately, 
the radial mode frequencies cannot distinguish between the 
two phases of stellar evolution - shell H-burning and core He- 
burning. This is hardly surprising, since the models inside the 
box on the HR diagram in either phase would have approxi- 
mately similar radii, and the large separation depends crucially 
on the radius of the star. However, the seismic information helps 
us to constrain the radius of the star to a much narrower range 
tha n that possible by the er rorbox on the HR diagram. 

iKallinger et alj (2008a) have carried out a seismic modelling 
study of € Oph also, but with important differences. Firstly, they 
have interpreted the observed peaks in the power spectrum of 
e Oph to be radial as well as nonradial modes. They identify the 
sharp narrow peaks in the power spectrum as long-lived nonra- 
dial modes, as compared to the broad Lorentzian envelopes, re- 
sulting in short-lived modes, w hich have been identified as radial 
modes by Bar ban et al.l (|2007). Secondly, they have matched the 
observed frequency values, and not the large separations, to the 
theoretical model frequencies. Lastly, they have used only shell 
H-burning models. Given these differences, it is not surprising 
that they obtained a slightly different result than ours. 

Little is known whether p-modes, radial or nonradial, can be 
excited in red giants to an amplitude high enough to be observ- 
ab le. A detailed t heore tical study of models of a UMa, observed 
by iBuzasi et all d2000l) with the WIRE sate llite, and of similar 
2M m odels on the lower giant branch, by iDziembowski et al. 
(2001) provides some insight on the oscillation characteristics 
of lower giant branch stars. Giant stars are characterised by an 
inner cavity that can support gravity waves (g-modes), and an 
outer cavity that supports acoustic waves (p-modes). Observable 
modes in the outer cavity are mixed modes, with a g-mode char- 
act er in the inner cavity and p -mode character in the outer cav- 
ity. Dziembowski et al. (2001) showed that such mixed p-modes 
can have substantial amplitudes, and that low degree modes with 
I = 2, 3, together with the radial /j-modes can be unstable. 
According to their models of lower giant branch stars, high am- 
plitudes in the outer cavity arise only for modes with 1-2. The 
excitation properties of /?-modes in more luminous red giants ly- 
ing around the middle of the red giant branch, which have very 
deep convection zones, such as e Oph, are mostly unexplored. 

A careful visua l examination of th e observed frequency 
spectrum of e Oph (Ba rban et alj 120071) reveals a clear comb- 
like structure with reasonably regular spacing of ~ 5.3 //Hz. In 
an adiabatic pulsation calculation of a theoretical model, a series 
of regularly spaced radial modes accompanied by a dense for- 
est of nonradial modes are obtaine d (see Fig. [8] ). The observed 
frequencies of e Oph, as given by iBarban et alj d2007l) . can be 
matched reasonably well to the radial modes. In principle, they 
can also be matched easily to many of the closely spaced non- 
radial modes as well, but the question remains as to why the 
majority of the nonradial modes are not observed. 
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Arbitrary radial order 

Fig. 8. Adiabatic pulsation frequencies of radial {€ = 0) and non- 
radial (I = 1,2,3) modes in a typical red giant model inside 
the e Oph errorbox on the HR diagram are shown. The frequen- 
cies of each degree are plotted in increasing order of magnitude 
against some arbitrary assignment of radial order (which can be 
simply treated as the serial number of the mode in the eigenspec- 
trum). The top row shows the frequencies of a 1.8M star inside 
the e Oph in the shell H-burning phase. The bottom row shows 
the frequencies of the same star when it again appears inside the 
e Oph errorbox in the later phase of stable core He-burning. The 
horizontal grey bands represent the observed MOST frequencies 
of e Oph with ±lcr error bars. 



A possible explanation of this may be provided in terms 
of the normalised mode i nertia E, defined according to 
IChristensen-Da lsgaard (2004j) as 



\ 2 dV 



M\£(R)\ 2 



(6) 



where £ is the displacement vector, and the integration is car- 
ried out over the volume V of the star. Most nonradial modes 
are strongly trapped in the core, with only a few modes being 
trapped in the envelope. The modes trapped in the core have high 
inertia and therefore unlikely to have enough surface amplitudes. 
On the other hand, the modes trapped in the envelope can have 
inertia as low as the radial modes, and therefore, be seen on the 
surface. Fig. |9] shows the normalised mode inertia for radial as 
well as nonradial modes in a typical model inside the e Oph er- 
rorbox. Only selected I - 1, 2, 3 modes have inertia comparable 
to the radial modes. However, a reliable estimate of the surface 
amplitude can only be obtained when nonadiabatic effects in the 
outer layers are taken into account. 
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Fig. 9. Normalised mode inertia of I — 0, 1 , 2, 3 modes as a func- 
tion of frequency is shown for a 1.8M Q star in shell hydrogen 
burning phase inside the e Oph errorbox. The individual eigen- 
frequencies are joined by lines for each degree (red filled circles 
with solid lines for £ — 0; green empty circles with dotted lines 
for I — 1 ; blue crosses with short dashed lines for £ — 2; ma- 
genta empty squares with long dashed lines for I — 3). 



In a recent theoretical work involving no nadiabatic treatment 
of the excitation mechanism, iDupret et al.l d2009l) have found 
that despite their low amplitudes, a selection of nonradial modes 
may still have appreciable heights in the power spectrum of red 
giants due to their long lifetimes, and thus be detected in obser- 
vations. However, the detection of such modes depends crucially 
on the evolutionary stage of the star on the red giant branch. 
Theoretical computations for intermediate red giant branch stars 
like € Oph predict much longer lifetimes for nonradial modes 
(k. 50 days) than radial modes (;£ 20 days) because of their larger 
inertia. But if the duration of observation is shorter than these 
lifetimes, as the case is for MOST observations of e Oph, it is 
not possible to resolve these nonradial modes in the power spec- 
trum. This implies that the modes attain smaller heights in the 
power spectrum and they become extremely difficult to detect. 

Further, in the higher part of the frequency domain, the de- 
tectable nonradial modes also appear with an asymptotic regular 
separation pattern similar to that found in main sequence stars. 
This means that the I — and i — 2 modes will appear close to 
each other, with the I = 1 modes occurring roughly midway be- 
tween them. For e Oph such a pattern implies that if we consider 
nonradial modes to be present in the spectrum, the large separa- 
tion (of i = modes, for example) would be almost double the 
value than that obtained by postulating only radial modes. Such a 
high value of the large separation (~ 1 1 //Hz) is inconsistent with 
th e position of e Oph on the HR diagram. However, according 
to lDupret et al.l d2009l) . the lifetimes of some 6 = 2 modes which 
are strongly trapped in the envelope are comparable to that of the 
radial ones. So these envelope-trapped 1-2 modes could indeed 
be detectable. Trapping of £ = 1 modes is less efficient, as found 
by Dziembowski et al. (2001) too, and hence they may not be 
observable. However, in the absence of any detailed calculations 
for the mode amplitudes and lifeti mes for the specific case of 
e Oph, in this work we have adopted Barban et al.l d2007l) 's inter- 
pretation of the frequencies as radial modes only. An alternative 
modelling analysis taking into account the possibility of nonra- 
dial mode s might lead to a differe nt set of model parameters, as 
found by Kallinger et alj d2008al) . for example. The theoretical 



w < 

X 



Shell H-burning o 
[Core He-burning 




j i i i i i i i i_ 



Fig. 10. The x\ v values are shown as functions of radius for all 
stellar models constructed within the e Oph errorbox. The red 
circles represent the shell H-burning and the blue crosses repre- 
sent the core He-burning models. The vertical dashed lines de- 
note the lcr interval for the interferometric radius. 



justification behind the presence of only a few specific nonra- 
dial modes among the possible dense spectrum of such modes 
requires further detailed study. 

Our interferometric measurements yield a value of the ra- 
dius of e Oph that is in close agreement with the radii of our 
best models obtained through seismic analysis. The interfero- 
metric radius was not used as a constraint to choose the best 
model, but was rather checked a posteriori against the seismic 
values. The range of possible values of radius of stellar models 
within the errorbox on the HR diagram for e Oph is 9R to 12/?©. 
But the seismic information, specifically the frequencies or the 
large separation, restricts the radius to a much smaller range. The 
best-fitting models using large separation comparison have radii 
of R » l0.5R e for both shell H-burning and core He-burning 
phases. This value is within 2cr of the interferometric radius. 
Even considering all models with ^ 2 < 1 constrains the range 
of radius to 10.4 < R/R e < 11.2 for frequency comparison, and 
to 10.2 < R/Rq < 11.6 for large separation comparison (see 
Fig. [Tol l. This range encompasses the much narrower limit for 
the radius set by the interferometric measurements. It is remark- 
able that despite the inherent uncertainties in the modelling of 
the outer layers of a star, the seismic analysis alone leads to a 
value of the radius that is in such good agreement with an in- 
dependent direct estim ate of the radius. The radius obtained by 
iKallinger etall d2008al) through frequency fitting, R = 10.8fl G , 
seems to be more removed from the interferometric radius than 
our seismic value is, although it is difficult to directly compare 
the two in view of the absence of error bars for the former. 

As mentioned above, in the present study the independent 
estimate of the radius was not used as an additional constraint 
for choosing the best seismic model. However, since the mod- 
elling is in no way influenced by the presence of the radius in- 
formation, we can check what would be the result if indeed the 
radius is used as a constraint. We use the 1 <x interval of the ra- 
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Table 4. Stellar parameters for the models with lowest j£ and xl values in either shell H-burning or core He-burning phase that 
have radii within +l<x of the interferometric radius are given. 



Minimised by 


Evolutionary phase 


Y 


Z a 


M/M e 


Age (Myr) 


log r eff 


log L/L 


R/Ro 


Ay 


Xs v 


XL OTXy 


Shell H-buming 


0.275 


0.012 1.8 


1.9 


985 


3.685 


1.732 


10.44 


3.38 


0.73 


Core He-burning 


0.270 


0.013 1.8 


1.8 


1322 


3.691 


1.754 


10.41 


3.42 


0.67 




3.700 3.695 3.690 3.685 3.680 
log T e „ (in K) 

Fig. 11. All the computed stellar models inside and around the 
e Oph errorbox on the HR diagram are shown. The red open 
circles represent the shell H-burning and the blue open squares 
represent the core He-burning models. The filled symbols rep- 
resent the models with x\ v ^ 0.8 for each phase of evolution. 
The big circle and the big square show the best models in the 
two phases. The solid line rectangle represents the adopted er- 
rorbox of e Oph. The dashed lines denote the loci of constant 
radii R = 10.32/? o and R = 10.46/?©, which are the + lcr bounds 
for the interferometric radius of e Oph. 



dius to restrict the position of the star on the HR diagram, along 
with the adopted values of luminosity and effective temperature. 
This leads to a trapezoidal area on the HR diagram (see Fig.Qj} 
as the errorbox for e Oph. The parameters of the models inside 
this smaller errorbox that fit the seismic data best are shown in 
Table [4] In this case, the minimisations according to frequen- 
cies and large separations yield the same best models in either 
phase of evolution. The minimum values of x\ Y are marginally 
higher than in the more general case (cf. Table [1]). However, the 
minimum xl values are significantly higher, indicating that al- 
though the observed large separations are quite well matched 
by these models, the absolute frequency values are somewhat 
shifted. This is wholly expected since the large separations are 
strongly influenced by the radius, even if the frequencies them- 
selves may be shifted due to inadequate modelling of the sur- 
face effects. In other words, restricting the radius value implies 
a strong constraint on the large separation, but not necessarily 
on the absolute frequencies. This is also borne out in the more 
general case (when the radius constraint is not used) in Fig.Qj] 
where the models with lowest x\ v values (< 0.8) in either phase 
of evolution lie in a broad band almost parallel to the interfer- 



ometric radius band, indicating a constant higher radius value 
(~ 10.55 + O.15/? ) common to all of them. The seismic values 
of the mass and radius of e Oph yield an average large separa- 
tio n value of 5.35 + 0.1 8 //Hz according to the scaling formula 
of Kieldsen & Bed dind (1 19951). which is completely consistent 
with the current data (iBarban et al.ll2007l) . 

It is also possible to determine the effective temperature of 
the star from its measured radius and photometric data. From 
Fig. Qj] itself, it is evident that the intersection of the luminosity 
limits and the interferometric radius ranges indicates a tempera- 
ture range of log T eS = 3.695 + 0.010 (4955 + 100K) which is 
consistent with our adopted range of effective temperature and 
has a similar uncertainty. Further, a fit of the available photo- 
metric data on e Oph (BVJHK bands) using tabulated Kurucz 
models yields log T eff = 3.691 + 0.002 (4912 + 25 K) for surface 
gravity values typical of stars in the relevant zone of the HR di- 
agram. Actually, a change of 0. 1 dex in log g makes a difference 
of only 1 K in T e g, while the uncertainty in the measured angu- 
lar diameter contributes about 4K in the error estimate. Again, 
these value of T e ff are complete l y cont ained in the range that we 
have used from de Ridd er et~ai] (|2006). Thus our adopted values 
of L and T e g are consistent with the independent measurements 
of parallax (Hipparcos) and the angular diameter (this paper). 

It is evident from this study that the seismic information 
alone can go a long way in constraining the most important stel- 
lar parameters of red giants. Even with a very limited data set, it 
was possible to obtain a reasonably narrow range of parameters 
for e Oph, and the radius estimate from the seismic modelling 
stands in close agreement with a completely independent in- 
terferometric measurement. However, the accuracy of the mod- 
els can be greatly enhanced by the additional information about 
the interferometric radius. An independent radius measurement, 
with a high precision such as provided by interferometry, helps 
in reducing the size of the errorbox on the HR diagram, mak- 
ing the task of searching for the best model easier. This is the 
first instance of the coming together of asteroseismology and in- 
terferometry for red giant stars, and clearly illustrates the huge 
potential of this combination in detailed studies of such stars. 
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